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Abstract. This paper is devoted to present an approximation of a Cauchy 
problem for Friedrichs’ systems under convex constraints. It is proved the 
strong convergence in Lj* of a parabolic-relaxed approximation towards the 
unique constrained solution. 


1. Introduction. The aim of this paper is to prove the convergence of a relaxation 
approximation of weak solutions to Friedrichs’ systems under convex constraints. 
The wcll-posedness has been established in [3] by means of a numerical scheme. 
We present here another way to construct such weak solutions thanks to a model 
that relaxes the constraints. We consider the following Cauchy problem: find W : 
[0, T] x 1” —> R m such that 

'd t W + £"= i BjdjW = 0 in ]0, T] x R n , 

< W(t, x) e K if (t, x) e [0, T] x R n , (1) 

W{Q,x) = W°{x) if x € R", 

where A is a fixed (i.e. independent of the time and space variables) non empty 
closed and convex subset of R m containing 0 in its interior, the matrices B :j are 
m x m symmetric matrices independent of time and space, and T > 0. The main 
difficulty is due to the constraints which introduce nonlinear effects to the linear 
Friedrichs’ system [5]. This type of hyperbolic problems has been introduced in [3] 
where a notion of weak solutions to problem 1 has been defined. 
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Definition 1.1. Let W° £ L 2 (W n ,K), and T > 0. A function W is a weak 
constrained solution of 1 if W £ A 2 ([0,T] x R”, K) satisfies 



(\W-K\ 2 d t cl ) + Y / (W 

j=i 


K\Bj{W 



dt dx 


+ ( \W°(x) — k\ 2 (I)(0,x) dx > 0, (2) 

J r 


for all k £ K and 4> £ C£°([0, T[xM n ) with (f>(t,x) > 0 for all (i,ir) £ [0, T[xR". 


We recall here the main result of [3]. 

Theorem 1.2. A ssume that W° £ L 2 (JBL n ,K). There exists a unique weak con¬ 
strained solution W £ L 2 ([0, T\ x R n , K) to 1 in the sense of Definition 1.1. In addi¬ 
tion, this solution belongs to C([0, T], L 2 (R", K)), and if further W° £ H 1 (W n ,K), 
then W £ L°° ([0, T ], H 1 (R n ,K)). 


As already mentioned, the well-posedness of problem 1 has been established in [3] 
thanks to a numerical method. We relax here the constraints W(t, x) £ K for a.e. 
(t,x) £]0, T[xR" as 


d t W e + Y,Bjd Xj W e 

3 =1 


P K {W e )-W e 

e 


( 3 ) 


where Pk denotes the orthogonal projection onto the closed convex set K, and 
e > 0 is a small parameter. Formally, if we multiply equation 3 by e, and let e tend 
to 0, we get that the “limit” of W e , denoted by W, satisfies Pk(W) = W, which 
ensures that W £ K. In addition, Definition 1.1 has been motivated in [3] by a 
formal derivation from the relaxation system 3. To see it, it suffices to take the 
scalar product of equation of (3) with W e — k, where k £ K is arbitrary. We then 
use the first order characterization of the projection which ensures that the right 
hand side is non-positive, to get the inequality of Definition 1.1. The purpose of 
this work is to rigorously justify these formal steps. 

The relaxation model presented here is very similar to viscous approximation of 
constrained models found in mechanics, and especially in plasticity. We start from 
the system of dynamical linear elasticity in three space dimension which can be 
written as 

f d t F + S7 x v = 0, ,... 

\ dtv + div cr = 0, 

for all t £ [0,T] and x £ R 3 . In the previous system, F(t,x ) is a 3 x 3 matrix 
which stands for the displacement gradient, v(t,x) £ M 3 is the velocity, i.e. the 
displacement time derivative, and a = fi (F + F T ) + A (tr F) I 3 is the symmetric 
Cauchy stress tensor (here A and /a are the Lame coefficients , and I 3 is the idendity 
matrix in R 3 ). This system can be rewritten (thanks to a change of variables - 
see [6]) in the Friedrichs’ framework as 


dtU + A\d Xl U + A2dx 2 U + A 3 d X3 U — 0, 

where U is a vector in R 9 (containing the three compononents of the velocity v, and 
the six components of the stress cr) and A 3 ,A 2 and A 3 are symmetric matrices. 

We now introduce the convex contraint coming from plasticity, see [11]. Indeed, 
the theory of perfect plasticity is characterized by the fact that the stress tensor a is 
constrained to stay inside a fixed closed convex set K of symmetric 3x3 matrices. 
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The total strain is then additively decomposed as the sum of (i) the elastic strain, 
denoted by e, which is still related to the stress by the linear relation a = Atre + 2/ie; 
(ii) and the plastic strain, denoted by p, whose rate is oriented in a normal direction 
to K at a. Summarizing, one has 

F + F t 

---= e +p, a = Atre + 2/ie £ K , d t p £ <9/r-(ct), (5) 

where dlx(.cr) denotes the subdifferential of Ik, the indicator function of K , at the 
point a. Using Fenchel-Moreau regularization of Ik (see [8]), the last condition in 
5 can be relaxed as 

d t p= - (P K (cr) - a), 
e 

where e > 0 is a viscosity parameter. We can now reformulate, at least formally, 
the dynamical problem of visco-plasticity (see [9, 11]) as 

[ dtU + A,d Xl U + A 2 d X2 U + A 3 d X3 U = P ^ U )~ U on ] 0 , T] x R 3 , 

4 U(t , x) £ K if (i, x) £ [0, T] x R 3 , ( 6 ) 

[l/(0,a:) = l7 o (a:) ifxeR 3 , 

where, again, U £ R 9 , and K = {u £ R 9 : a £ K }, and Ai, A 2 and A 3 are the 
same matrices than in the elasto-dynamic case. As e tends to zero, one expects the 
solution to 6 to converge to that of the model of perfect plasticity (see [12] in the 
quasistatic case). 

Notation. In the sequel, we denote by (|) the scalar product of L 2 (R n ,R m ) and 
by (;) the canonical scalar product of R m (and |.| the associated norm). Also, to 
shorten notation, we write L 2 x (resp. H^ x ) instead of L 2 (0, T; L 2 (R”,R m )) (resp. 
.ff 1 ((0,T) x R",R m )), L 2 instead of L 2 { R n ,R m ). 

This paper is organized as follows. In the first section, we establish the existence 
and uniqueness of the relaxation model thanks to a parabolic approximation. In the 
second section, we prove that the relaxed solution W e to 3 satisfies the inequalities 
of Definition 1.1. Finally, to get the existence of a solution as a limit when e tends 
to zero of relaxed solutions W e , we prove the strong convergence of the sequence 
(bF e ) e >o in the space L 2 ((0,T) x w), where w is a open bounded subset of R", to a 
weak solution to the contrained Friedrichs’ systems. 

2. Parabolic approximation. In order to find a solution to the relaxation prob¬ 
lem 3, we use a parabolic type regularization. To this aim, we consider a classical 
sequence of modifiers in R”, denoted by (p^)r;>o- 

Theorem 2.1. Let W° £ Id 1 (R", K). For every e > 0 and r] > 0, the system 

f d t W £ , v - r/AW^ + BjdjWe,,, = Pk{W ^~ W ^\ 0 n ]0, T) xR", 

| 1=1 ^ ' 
[w ej7) (0,x) = W°(x) * p v , ifxe R n , 

admits a unique solution We )T? with the following properties: 

W £tV £ L 2 (0,T; 7J 2 (R", R m )), d t W^ £ L 2 (0, T; IJ^R", R m )), 

and 

d tt W etr ,eL 2 (0,T-H- 1 {M n ,W n )). 
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Furthermore, we have the following estimates 


sup \\W e , v (t )\\ 2 L 2 
0<t<T x 

< 

II rr II LI ’ 

( 8 ) 

sup \\W e , v (t)W Hl 

0 <t<T 

< 

C € \\W°\\ Hl , 

(9) 

sup \\d t W et r,{t )\\ L 2 

0 <t<T 

< 

c t \\w°\\ Hl , 

(10) 


for some constant C e > 0 independent ofr). 


Proof. The proof essentially follows that of Theorem 1, Part II, Section 7.3.2 in [4]. 
Let X = L°°(0,T-, H 1 (W 1 and V £ X. We consider the problem 


jdtU- r)AU = Pk ^ —— -^TBjdjV, on ]0, T\ x 1R”, 
\u(0,x) =W°(x) * p v , if x £ R n . 


( 11 ) 


Since 0 £ K, we have the following inequality 


Vfc £ R m , \{P K {k)-k)\<\k\ 

which shows that BjdjV £ L 2 ( 0, T; L 2 (R", R m )). Using the theory 

of parabolic equations, we get that equation 11 admits a unique solution U with 
U £ L 2 (0,T;H 2 (R n ,R m )) and d t U £ L 2 (0, T; L 2 (R", M m )). Let V £ X and U be 
the solution to 11 associated with V. The function U = U — U is a solution to 


d t U — yAU = 

< 

17(0, x) = 0, 


Pk(V)-Pk(V)-(V-V) 

€ 


B j9j(V-V), on ]0, T\ x 


( 12 ) 

Thanks to the theory of parabolic equations, we have the following estimate, we 
obtain the following estimate 


ess-sup 
o <t<T 


m 

<C(v) 


HI 


Pk(V)-P k {V)-(V-V) 


- £ B A(V - V) 

3 =1 


< C(rf) max ( -, ||B 


V-V 


Li hi 


< C{rj) max ( ||5j|| ) VTess-sup V(t) — V(t) 

\ e ) 0 <t<T 


HI 


Therefore, the mapping 


ip : 


jx -A {U £ L 2 H 2 and d t U £ L 2 t x ) C X 


V^U 


is Lipschitz-continuous with Lipschitz constant bounded by C(rj) max \\Bj ||) y/T. 

We now divide [0,T] into sub-intervals [0, Ti], [Ti, 2Ti], [2Ti, 3Ti],..., [TVTi, T\ 
such that 

C{rf) max ||Bj||^ ^/max(Ti, T - IVTi) < 1. 
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For every i £ {0,... AT — 1}, the Banach fixed-point Theorem ensures the existence 
and uniqueness of a solution to 7 on the interval [iT\, (i+l)Ti] (with initial condition 
W° * p v if i = 0, and W (:V (i r I\) if i > 1 obtained at the previous step). We 
then obtain a solution W e>ri on the entire interval [0,T] by gluing the solutions on 
each sub-intervals, so that W e ^ £ L 2 (0,T; ff 2 (R n ,R m )). According to the initial 
condition on each sub-intervals, the function t H > VF £j 7 ) (t) is continuous in L 2 at 
every t = iT\, so that dtW e , v £ L 2 (0,T;L 2 (R n ,R rn )), and, in particular W e , v £ 
H 1 (]0,T[xR n ,R m ). 

To obtain the announced regularity, we use the following result (whose proof 
relies on the chain rule in Sobolev spaces). 

Lemma 2.2. Let U £ i7 1 (]0,T[xR n ,R m ). Then the function Pk ° U — U also 
belongs to i7 1 (]0, TfxK",® 1 ™), and there exists a constant M > 0, independent of 
U, such that 

\\PkoU-u\\ h1x <m\\u\\ h1x . 

This result ensures that PKl - n ‘-P~ u e - T| £ i7 1 (]0,T[xR”,R m ) and then, using 
again to the regularity theory of parabolic equations, we obtain that 

W 6 ,„ £ L 2 (0, T ; H 2 (R n 1 R m )), d t W e , v £ L 2 (0, T; H\R n , R m )), 

and 

d t tW e ,n £ L 2 (0,T-,H- 1 {R n ,R m )). 

Now we derive the estimates. We are going to use the following result (see [4]). 

Lemma 2.3. Let U £ L 2 (0,T,H 1 (R n ,R m )) with d t U £ L 2 (0, T; H~ 1 (R n , R m )). 
Then, the function 

t \\U(t)\\ Ll , 

is absolutely continuous, and for a.e. t £ [0,T], 

1 QtoIIlj) =(u(t),d t u(t)) H ^. 

Applying this result to W eiV , we get that for a.e. t £ [0,T], 


d 

dt 


0 II^C^WIIl 2 ) — {W e> rj(t)\dtW e>r j(t)) L 3 


= w eiV (t) 


^w - £ B A W, M + ^ W , M \, (13) 


1=1 


where we used the fact that W e , v is a solution to the partial differential equation 7. 
Since 0 £ K, we have 


(w etri (t)\p K (w^m-w e , v (t)) L2x 

= {p K {w e , v m i p K (w e , v m - w e , v (t)) - ii p K (w e , v m - w^ml* < «• ( i4 ) 

On the other hand, if v £ C(N(R",R. m ), an integration by parts shows that 

(v | fjA«) L2 = -77 \\Dvf L 2 < 0, 

and 

n 

B J d J v 

f=i 



V 
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since the matrices Bj are symmetric and independent of the space variables. By 
approximation, these formulas are true for v £ H 1 (as well, and in partic¬ 
ular, 


W e , v (t) 


E 

3=1 


BjdjW tiV (t) + rjAW etri (t) } <0. 


(15) 


Gathering 13, 14 and 15, we obtain that 


d 

dt 


\W e J 


ll 


< 0 , 


and using Gronwall’s Lemma, we derive the first estimate 8. 


sup ||W CjT7 (t)|| 
0 <t<T 


2 

L l 


w°*p4 2 li < \\w° 


2 

Ll ’ 


(16) 


We apply the same argument to the spatial weak derivates d k W f . v of W £tV . De¬ 
riving the partial differential equation 7 in the sense of distribution, we infer that 

d t d k W e>v - V A d k W^ + J2 BjdjdkW. r e , v = dk PK{W ^)- w ^ K 

3=1 

The previous equality actually holds in L 2 (0, T, L 2 (M. n , R m )) thanks to the regular¬ 
ity of W etV , and we can apply Lemma 2.3 to obtain that Efor a.e. t £ [0,T], 


j t (\ \\d k W^(t)\\ 2 Li ) = (d k W e , v (t),d t d k W e , v (t)) Li 


= ( d k W e , v (t) 


0 PK(W e ,r,)(t) - We, v (t ) ^ ^ , A n „ r E 

d k -2^ Bjdjd k W e ^{t) + r]Ad k We,r,(t) 


3=1 


Arguing as in 15, we get 


lll^j < ^W^t) 


P K (W e<ri (t))-W t>v (t ) 


<-\\d k W^{t)\\l^ 

6 x 

where we used Lemma 2.2. Using again Gronwall’s Lemma, it yields 

sup \\d k W e ^(t )\\ 2 L 2 < exp \\d k W° * prj ||E 

o <t<T x \ e / 


< exp 


^2TM ^ 


I d k W° 


I Ll ’ 


(17) 


which completes the proof of estimate 9. 

We finally derive the last estimate for d t W f , rr Again, we derive the partial 
differential equation 7 with respect to t in the distributional sense to get 

d tt W.^ - vAd t W e , v + BjdjdtW^ = d t PK ^ W ^ ~ 
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in L 2 (0,T,H 1 (R”,R m )). Using again Lemma 2.3, we get that for a.e. t £ [0, T], 
d (\ 


dt V 2 


||dtW ej!7 (t)|| i 2 ) — (dtWz^it) | dttW e<ri (t)) H1 


HZ 


= ( m^t) 


AO IT/ CA D a O w CA o -Pff(U / e,rj)(i) — We,?j(i) 

rjAd t We t ri(t) / , Bjdjd t W e}71 (t) -{- Of 


7=1 


As before, we infer that 




y: Bjdjd t W eiV (t) \ =0, (5 t We,^(t) I r]Ad t W et ri(t)) H 1 >ff -i < 0, 


7=1 


Id 


which shows, thanks to Lemma 2.2, that 


d (l 


J t ( 2 WdtW^mli ) < 7 \\d t W^(t)\\ Ll II d t {P K {W^){t) - W e , v (t))\\ L 2 


<^\\d t w e , v m 2 Ll - 


At this point, we would like to use Gronwall’s Lemma. To do that, we need to 
know the value of d t W ejV a.t t = 0. To this aim, let us take a test function cj> £ 
C£°(M x R n , R m ) with <j)(T, •) = 0. On the one hand, according to Fubini’s Theorem 
and Green’s formula on [0,T], Ewe have 


/ / (d t w e , v -M) = 

10 J 1" J K" 


- <^^,^(0,-);0(o,-)> - [ (dttW e , v -A) 

Jo 


since d t W e>7J £ H 1 (0,T;H 1 (R”,R m )) and </> is smooth. On the other hand, ac¬ 
cording to equation 7 


io J R" 


(dtW^dtcfi) 


f ( (rjAW^-j^BjdjW^ 
Jo JWL n \ j=1 


Pk(W^)-W £ , v . 


\d t 4> 


r?AW e ,„ - ^ B A W e,r, 
7=1 


PKiW^n) - W e , v 


(’ 


- ( 7]AW° B.djW 0 * p v + Pk ^ W °*P^ ty °*^ ;0(O, .) 


7=1 


Since </>(0, •) is arbitrary, we obtain that 

d t W e , v (0, •) = r)AW° * Pr ,-J2 B A W ° * Pv + Pk ^ W ° * Pn l ~ W °* P \ 

7=1 6 

We are now in position to apply Gronwall’s Lemma which implies that 


( 2 TM 

— 

e 


||9tWe,^(0, *) llx,2 ■ 
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Since, 


c 


A w°*p v \\ 2 < - ||VW' 


V 


I LI ’ 


for some constant C > 0 independent of 77, we deduce that 


sup \\dtW e ,r,(t)\\ L2 <C e \\W° 


0<t<T 


I Hi ’ 


where C e > 0 is another constant independent of 77, which completes the proof of 
the last estimate 10. □ 


3. Approximation of the convex constraints. We now consider the relaxation 
problem 

d t W ( + V BjdjWt = Pk ^ W ^ ~ We , on ]0, T] x R", 

U e ( 18 ) 

W e (f),x) = W°(x), if x £ R". 

Thanks to Theorem 2.1, we will construct the solution to the previous problem as 
the limit of the solution of the parabolic problem 7 when 77 tends to zero. 

Theorem 3.1. There exists a unique solution W e £ .H^QO, T[xR", R m )) to 18 
satisfying, for all <j> £ C£°(R x R”,R m ), 

f f / d t W e + BjdjW e \<j\ dx dt = [ [ ( Pk ^ ~ We ■ A dxdt, (19) 

Jo J R n \ - =1 j Jo J R" \ e / 

and W e ( 0,-) = Wq in L 2 (W n ,W m ). In addition, 

sup \\WM\l* < \\W °\\ 2 la • ( 20 ) 

0 <t<T 

Proof. Thanks to Theorem 2.1, the sequence (W e ^) T) >0 is bounded in the space 
17 1 (]0,T[xR”,M m ). We can thus extract a subsequence (not relabeled) such that 

W 6 ,„ ->• W e weakly in T[xR n , R m ). 

In particular, 20 is a consequence of 8 by the lower semicontinuity of the norm 
with respect to weak convergence. Since the embedding of il 1 (]0, T[xM", R. m ) into 
Lf oc ([0,T] x M",!!" 1 ) is compact (cf [1]), we deduce that 

W e>v W e strongly in L 2 (]0,T[x] - R, R[ n , R m ), 

for each R > 0. Let <f> £ C£°(R x R”,R m ) and R such that the support of is 
contained in [-R, R] n+1 . Since W tyV is a weak solution of 11, we have 


n 

to JR n 


n n 

( {d t We, v ; 4>) + 77 ^ (9jW e ^; dj4>) + ^2 (BjdjW'rf <f>) ) dx 


dt 


3 =1 


3 =1 

r T f / P K {w e , v )-w e , v . 


to J R n 


v 


\(j)) dx dt. 


Using the weak convergence, we infer that 

rT 


'0 7 R" 


3 =1 


({dtW €>r ,;<f>} +r]'^2 ( d jW e ,n',dj4 >) + y^ t {BjdjW^r,‘A) ) dxdt 

' ' j=i 

n 

( ( dtW e -,(j)) + ^ (BjdjW e -, </>)) da; dt. 


n 

to JR n 


3 =1 
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On the other hand, the strong convergence yields 
^ P K (W e , v )-W e>ri ^ ^ /\ P K (w t )-w, 

and consequently, we obtain that 

r T 


;</>) strongly in L 1 (]0,T[xR"), 


/o J R n 


d t W e + ^2 BjdjW e -, (pj dxdt = J (EYYA — ^--^jdxdt. (21) 


We next focus on the initial condition. We take </> £ C“(] — oo,T[xI 
particular <fi(T) = 0). An integration by parts shows that 

— [ f (W £jr) ; dx dt + f (W° * p v (x)\ </>(0, x)) dx 


m TO)ra 


) (i 


(9(144^; (j>) dxdt. 


/ o Jr™ 


Letting r] tend to zero, and using 21 leads to 

rT r n 

( - (W4; d t (f>) + Y (BjdjW^c/)) ) dx dt 


/0 JR™ 


1=1 


+ 


{W°(x);<l>{0,x))dx= f f 
Jo Jr 


^ P K (W e )-W e 


\(j)\dxdt, (22) 


since W° * p v —> W° strongly in Lj J oc / TD,n 


r 2 ron^ ]R m ). We now integrate by parts in 21, 
using the fact that W e £ Lf 1 (0,T, L 2 (R”,M m )), 

r T 


/ 0 JE“ 


( <W 6 ; d t 4>) + Y { B AWe; <t>) ) dx dt 


t=i 


+ [ (W e (0, x)\ 0(0, x)) dx = f f / k(W c ) — dxdt. (23) 
Jr™ Jo Jr™ \ e / 

Using the equations 22 and 23, it gives that 

f (W e (0, x); 0(0, x)) dx = f (W°(x)',(j)(0,x)')dx, 

Jr™ Jr™ 

and consequently the initial condition is satisfied in L 2 (R”,R m ). 

It remains to show the uniqueness. Let us consider two solutions W e and W e as¬ 
sociated with the same initial condition W°. Using the partial differential equation 
18, we obtain that 


n 

dt (w e -We)+Y B o d J ( M b - W 6 ) = 


P K {W t ) - W £ - P K (We) + W e 


3=1 


As already observed, we know that for a.e. t £ [0,T], 

r. n 

/ Y ( B 3 d 3 w &) - ^(*); w e (t) - w e (t)) dx = o, 

./lUn Z - \ / 


and also 


1 d_ 

2 dt 


3=1 


W e (t) - We(t) L2 = ( d t W e (t ) - d t W e (t) W e (t) - We(t) ) . 
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Consequently, we have that 


1 d 

2 dt 


W e (t)-W e (t) 


Ll 


^ P K {W e )(t) - W e (t) - P K (W e )(t) + W e {t) 


w e (t) - W e (t) 
2 


< - 

e 


W e (t)-W e (t) 


2 

Ll 


since the projection is 1-Lipschitz. Gronwall’s Lemma thus implies that W £ = W f 
since they satisfy the same initial condition. As a consequence of the uniqueness, 
we deduce that there is no need to extract a subsequence from (W £tV ) v> o to get the 
convergences as g —> 0. □ 


4. Convergence of the relaxed formulation. In this section, we first show that 
the solution W e to the relaxation problem 18 satisfies the inequality of Definition 1.1, 
and then we prove that we can pass to the limit in this inequality to get a solution 
to the initial problem 1. 

Lemma 4.1. Let W, be the unique solution to 18. For all k £ K and for all 
(j) g 1 / L 1 ’°°((—oo, T) x R™,R + ) with compact support in } — oo,T[xR n , one has 



(iWc-rfdtfi+Y.iWe 

3=1 


k; Bj (W e 



dt dx 


+ ( \W°(x) — k| 2 ^(0, x) dx > 0. (24) 

J R n 


Proof. Since W e is a solution to 18, we know that 



d t W f + Y, BjdjWe 


3= 1 


Pk(W.)-W., W '_ \ 


dx dt = 0. 


By the first order characterization of the projection, one has 


/ Pjc(wg~ w e 



( t , x) < 0 


for a.e. (t,x) £ [0,T] x M n . On the other hand, since W e £ H 1 (0,T;L 2 (W l ,K m )) 
and W e £ L 2 (0,T;lL 1 (R”,K m ))) we can integrate by parts to obtain the desired 
result. □ 


Remark 1. Let us stress that, although the function W e satisfies the same inequal¬ 
ity than the weak constrained solution, it is not a weak constrained solution in the 
sense of Definition 1.1 because it does not a priori belong to K. 

To get a weak constrained solution from the sequence of solutions (W e )e>o to 
the relaxation problem 18, we need to pass to the limit as e —> 0 in the previous 
inequality. This is the purpose of the following result. 

Theorem 4.2. For every bounded open set to C R”, the sequence ( W e ) e> o converges 
strongly in L 2 ((0,T) x w,R m ) to some function W which is a weak constrained 
solution to problem 1. 
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Proof. Let uj be an open bounded subset of R". We are going to prove the existence 
of a subsequence of (W e ) e> o (associated with the same initial data W°) which 
converges in L 2 (0, T, L 2 (w, R m )). We will use the following compactness criterion 
(see [ 10 ]). 

Theorem 4.3. Let B be a Banach space. A subset F of L 2 (0,T\ B) is relatively 
compact if and only if both conditions are fulfilled: 

• the set | J ( * 2 f(t) dt: f £ F^ is relatively compact in B for all 0 < t\ < < T; 

• we have 

sup \\r h f - f \\ L 2 (0 T _ h B) -> 0 , 
where r h f : (t, x) H- r h f{t, x) := f(t + h, x). 


We are going to apply this result to F = {W e \J)e where WW, is the restriction 
to [0, T] x u of W e . We first show that the set F = j J t * 2 W e ^(t, •) dt : e > 0 j is 

relatively compact in L 2 (oj, R m ) for all 0 < t\ < t 2 < T. We hrst observe that 
F is bounded in L 2 (w,R m ) by (t 2 — ti) ||W °|| i2 thanks to estimate 20. To show 
that F is relatively compact in L 2 (w,R m ), it is enough to check the validity the 
Riesz-Frechet-Kolmogorov compactness criterion (see [2] remark 13 page 74), i.e., 


lim sup 

h ^-o g>o 


rt 2 


(W e (t,x + h) 


W e (t , x)) dt 



(25) 


Note that W t ,h ■= W e (-, ■ + h) is a solution to the problem 18 associated with the 
initial condition W°(- + h). Consequently, since the projection is 1-Lipschitz, we 
have for all </> £ W 1 , 00 ((—oo, T) x R”,R + ) with compact support in ] — oo,T[xR ra , 



(dt(W e , h - W e ) + '52B j d j (W e , h - W e ),W e , h - W^jpdtdx 

= - f [ {P K (W e , h ) - P K (W e ) - (W e , h - W e ), W e , h - W € )<pdt dx < 0, 
e Jo JR" 


which implies that 



n 

(| W e , h - W £ \ 2 d t y> + (We,h - W e ,Bi(W e>h - W 6 )> d Xi v) dtdx 

i— 1 

+ [ \W° - W%\ 2 (x)cp(0,x)dx > 0, 
J i" 


where = W°(- + h). Let r > 0, we define the function p as 


f + V nLT ’ if t € [0,T] and r < |x| < r + nL(T — t), 
<p(t,x) = < if £ £ [0,T] and x £ 73 ( 0 , 7 -), 

\ 0 , otherwise, 


(26) 


where L is the maximum of the spectral radii of the matrices Bi. We claim that for 
a.e. (t, x) £ [0, T] x R ra 

(\W t ,h - W e \ 2 d t ip + J2{We, h - We, Bi(W e ,h - We))d Xi ^j (t,x) < 0. 
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This inequality is obviously satisfied as soon as x £ 5(0, r) and t € [0, T]. In the 
case where r <\x\< r + nL(T — f), we get for all 1 < i < n 

(W e , h - W e ,Bi(W e , h - W e )) (t,x) > -L\W e , h - W e \ 2 (t,x). 

Consequently multiplying by d Xi ip, it yields 


\^2^ We ' h ~ We, B i(W eth - W e )) d Xi ipj (t,x) <-\W eth -W € \ 2 (t,x)d t ip(t,x), 

and the inequality is true also in that case. According to 26, we obtain that 

f \W° -W°\ 2 (xM0,x)dx 
J R n 

rT /* n 

>-/ / {\W e , h -W t \ 2 d t ip + Y^(W e>h -W e ,Bi(W e , h -W e )}d Xi <p)dtdx. 

JO JB{ 0,r) ^ 

Thanks to the definition of tp, we get 

f f \W e , h - W e \ 2 dx dt < T [ \W° — W°| 2 dx, 

Jo J B(0,r) J B(0,r-\-nLT) 

and the regularity of Wo together with [2, Proposition 9.3] yields 

f [ \W e:h - W e \ 2 dx dt < T\h\ 2 [ \VW°\ 2 dx. 

Jo Jb( 0,r) ’ J R" 

Therefore, 25 holds, and consequently the set T is relatively compact in T 2 (w,R m ) 
for all 0 < t\ < t 2 < T. 

It remains to show that 

lim sup \\r h W e — W e \\ L2 ^ T -h-L 2 (ui,R m )) 0- 

h ^° e>0 v > V • n 

For all (j) € W 1,00 (] — oo, T — /i[xR", K + ) with compact support in ] — oo, T — /i[xR”, 
one has 


rT-h 


/o JR n 


= - [ [ ( P K (r h W e ) - P K (W e ) - (T h W e - W e ),r h W e - W e )<pdtdx < 0, 

e Jo J 8” 

since the projection is 1-Lipschitz. Arguing as before, we obtain 


{dt(r h W 6 - W € ), ( r h W e - We))ipdtdx 


rT—h 


J (y.BjdjfaW' - W £ ), ( T h W e - W e )^j tpdtdx 


rT-h 


| ThW e — W e \ 2 dtipdtdx 

rT-h r n 


rT-h r n 

+ / / Y. ( T h W e - W e ,Bi(r h W e -W e ))d Xi ipdtdx 

Jo JR " i=1 

+ f \W e (h,x) — W o (x)\ 2 tp(0,x) dx > 0. 
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Using a similar test function 

T T-~h + n L(T-h) ’ ^ t £ [0, T — ft] and r < |x| < r + nL{T - t), 
if t £ [0, T — ft] and x £ -B(0, r), 


‘KM) = { 


0 , otherwise, 

we get that for a.e. (t, x) £ [0, T — ft] xl", 


KW e - We| 2 d t </> + ^ (r^W, - W e , Bi{r h W e - W e )> d Xi y (t, x) < 0, 


i= 1 


and then 

rT—h 


[ [ \r h W e - W e \ 2 dx dt 

Jo Jb( 0,r) 


<(T-h) [ 
J E 


|W e (ft, x) — W°(x)\ 2 ip(0, x) dx 

B(0,r+nL(T-h)) 

<T [ | fU e (ft, x) — W°(x)\ 2 dx. (27) 

J B(0,r-\-nLT ) 

The conclusion then follows from the following result whose proof is very close to 
that of [3, Proposition 7}. 

Lemma 4.4. For all £ £ C£°(R",R + ), one has 


lim sup [ \W e (h,x) — W°(x)\ 2 ((x) dx = 0. 

h ^-o e>0 J R« 


According to Theorem 4.3 and estimate 20, the sequence (W e )e>o admits a sub¬ 
sequence (not relabeled) which converges strongly in L 2 ([0, T]\L 2 (ui, M m )) to some 
W and weakly in L 2 ([0, T}-L 2 (R n , R m )) to some W £ L 2 ([0, T]; L 2 (R", R m )). By 
uniqueness of the limit, we infer that W = W £ L 2 ([0, T]; L 2 (R ra , R m ). Let us take 
a test function <f> £ C(?°(]0, T[xR") in 19. Multiplying this inequality by e and pass¬ 
ing to the limit as e — > 0 yields W = Pk(W) a.e. in ]0,T[xR n which shows that 
W £ L 2 ([0,T]; L 2 (R”, A")). Finally, passing to the limit as e —> 0 in 24 shows that 
W is a solution in the sense of Definition 1.1 to the problem 1. Note finally that, 
by uniqueness of the solution to 1 (see [3, Lemma 9]), there is no need to extract a 
subsequence to get the above convergences as e —> 0 . □ 


The construction of the solution W to 1 rests on the assumption that the initial 
data Wo £ H 1 (R", K). Let us now explain how to construct a solution W when 
Wo only belongs to L 2 (R”, K). We use here the following result whose proof can 
be found in [3] 


Theorem 4.5. Let W° and W° £ H 1 { R",A"). We denote by W (resp. W) the 
solution in L 2 (0, T; L 2 (R ra , A")) to problem 1 in the sense of Definition 1.1 associated 
with W° (resp. W°). Then, W and W belong to C([0, T]; L 2 (R”, K)), and, in 
addition, we have the following estimate 


Vi £ [0,T],Vr > 0, W{t, ■) - W{t, ■) 


•) 

< 

O 

o 


L 2 (B(0,r)) 



L 2 (B(0,r+nLT)) 


where L is the maximum of the spectral radii of the matrices Bi. 
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By mollification, let us construct a sequence (WjJtgN such that W% £ H 1 (M n , K) 
for all k £ N, which converges to W° in L 2 (R n , K). The estimates of Theorem 4.5 
imply that 

sup ||W fc (t, -) - Wi(t ,-)|| ia(R ») < \\W° - W?\\ , (28) 

te[o,T] 1 ' 

where Wk (resp. Wj) is the solution to 1 associated with the initial condition 
Wj? (resp. Wi). It follows that the sequence {Wk)ke n is of Cauchy type in 
L°°(0,T; L 2 (R",R m )), and therefore it converges strongly in L°°(0,T; L 2 (R ra ,R m )) 
to some function W £ L°°(0, T; L 2 (R ra , R m )). Thanks to the strong convergence, 
we find that W satisfies the inequality 2. In addition, since Wk = Pi<(Wk) for all 
k £ N, we deduce that W = Pk(W) which ensures that W £ L°°(0, T; L 2 (R n , K )). 
The following result has thus been established. 

Theorem 4.6. Let W° £ L 2 (R n ,K), then there exists a unique solution W £ 
L°°(0, T, L 2 (R ra , K)) to 1 in the sense of Definition 1.1. 

5. Conclusion. In definitive, the relaxed problem 18, that was used in [3] to derive 
formally a definition of weak solutions of hyperbolic constrained problems, is in fact 
a rigorous way to construct weak solutions of hyperbolic constrained problems. It is 
worth noting that this relaxation procedure is deeply related to viscoplastic models. 
In order to fully apply this theory to mechanical problems, one should consider 
problems that are posed in bounded spatial domain. To do so, a new formulation of 
weak solutions to Friedrichs’ systems posed in bounded domains is proposed in [7], 
without constraints. It remains now to investigate the interactions between the 
boundary conditions which are considered in [7] and the convex constraints. 

Acknowledgments. The authors wish to express their thanks to Bruno Despres 
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mention that they have benefited from the unrivaled working atmosphere of the 
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